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We calculate the self-energies and the spectral densities of longitudinal and transverse gluons at 
zero temperature in color-superconducting quark matter in the color-fiavor-locked (CFL) phase. We 
find a collective excitation, a plasmon, at energies smaller than two times the gap parameter and 
momenta smaller than about eight times the gap. The dispersion relation of this mode exhibits a 
minimum at some nonzero value of momentum, indicating a van Hove singularity. 

■ PACS numbers: 12.38.Mh,24.85.-(-p 

o ■ 

, At asymptotically large quark chemical potentials and sufficiently small temperatures T, quark matter is a color 
superconductor jj. While there are, in principle, many different color-superconducting phases, corresponding to 
^ [ the different possibilities to form quark Cooper pairs with definite color, flavor, and spin quantum numbers, for 
—5 • three quark flavors and quark chemical potentials much larger than the strange quark mass, the ground state of 
color-superconducting quark matter is the so-called color-flavor-locked (CFL) phase 

In the CFL phase, the SU{3)c x SU{3)f color and (vector) flavor symmetry of QCD is broken to the diagonal 
CN ' subgroup SU{3)c+f. Consequently, all eight gluons become massive due to the Anderson-Higgs mechanism. The 
J> , situation is quite similar to that in the 2SC case, although there SU{3)c is broken to SU{2)c, and only five gluons 

■ become massive, while the three gluons of the residual SU{2)c color symmetry remain massless. These general 
5l5 [ symmetry considerations can be confirmed by an explicit calculation of the gluon Meissner masses in a given color- 

' superconducting phase. The Meissner mass (squared) of a gluon with adjoint color a is the zero-momentum limit of 
' the transverse component of the gluon self-energy at zero energy, limp^^o ^aaiPo — 0, p). 
\^ , At asymptotically large /z, the QCD coupling constant g 1, and the gluon self-energy is dominated by the 
' contributions from one quark and one gluon loop. The quark loop is ~ g^^J^^ while the gluon loops are ^ g^T^- 
,, Since the color-superconducting gap parameter is ~ ^exp(— l/gr) <SC 0|, and since the transition temperature to 
O • the normal conducting phase is (f>, for temperatures where quark matter is in the color-superconducting phase, 
T ^ Tc <^ H, the gluon loop contribution can be neglected. Following this line of arguments, the gluon Meissner masses 
^ have been computed for the 2SC phase in Ref. Q and for the CFL phase in Ref. |^ 0] . The full energy- momentum 
(— I [ dependence of the one- loop gluon self-energy has also been computed, but so far only for the 2SC phase The 
corresponding calculation for the CFL phase is the goal of the present paper. 

For the 2SC phase, the derivation of the gluon self-energies and propagators was explicitly discussed in Sec. II of 
Ref. 0, and we briefly remind the reader of the main steps. Just like in other systems (ordinary superconductors, 
\^ ' the standard model of electroweak interactions) where gauge symmetries are broken, there are Nambu-Goldstone 
5^ [ fluctuations which mix with the longitudinal components of the gauge field. For a particular choice of gauge ('t Hooft 
gauge), the Nambu-Goldstone modes decouple from the longitudinal components of the gauge field. Furthermore, for 
a particular choice of the 't Hooft gauge parameter, the Nambu-Goldstone modes can be eliminated and the gauge 
field propagator is explicitly 4-transverse. 

In our case, i.e., for the CFL phase, the computational steps leading to the gluon self-energies 11^^ and the gluon 
propagators A'^^ are quite similar to those in the 2SC phase 0, if not even simpler because all eight gluons are 
affected similarly by the breaking of the gauge symmetry, 

Il^^-^{P)^Sa,U^'^{P), AZ{P)^5abA^^''{P)., a,6=l,...,8. (1) 

Therefore, we do not give the details of this lengthy but straightforward calculation and just quote the final result. 
The final answer for the propagator of transverse gluons reads [cf. Eq. (56) of Ref. @] 
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FIG. 1: The imaginary parts of (a) II'''^, (b) —if^pi, (c) H^, (d) H*, and (e) n"" as a function of energy po for a gluon momentum 
p = 4(p. The sohd hnes are for the CFL phase, the dotted lines correspond to the HDL self-energy. 



while for longitudinal gluons we have [cf. Eq. (57) of Ref. 

A"°(P) = 1 . (3) 

p2_no"(P) 

Here, the transverse and longitudinal components of the gluon self-energy are computed from projections of H^'^ , 

n*(p) = ^{s^^ -py)n^^{p), (4a) 

-on. ^ 9 no"(p)tf(p) - [no^(p)pJ' , ^ 

where 

n'^{P)^p'Il'^P)p^ . (4c) 

The computation of the individual components and projections H™, 11°*, 11*, and 11^ is lengthy and thus deferred 
to the appendix. Note that the particular form of the longitudinal self-energy fl'^'^ arises from decoupling spatially 
longitudinal and time-like gluon degrees of freedom, see Ref. @]. 

In Fig. ^ we show the imaginary part of several components of the gluon self-energy for a gluon momentum p — Acj) 
as a function of the gluon energy po- For comparison, we also show the corresponding results for the gluon self-energy 
in the "hard-dense loop" (HDL) limit, Hg"^, cf. Eqs. (63a), (65a), (69a) and (69b) of Ref. The imaginary parts are 
quite similar to those of the 2SC case, cf. Fig. 1 of Ref. Q. There are subtle differences, though, due to appearance 
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of two kinds of gapped quark excitations, one so-called singlet excitation with a gap 0i, and eight so-called octet 
excitations with a gap = ^ In weak coupling, the singlet gap is approximately twice as large as the octet gap, 
01 ~ 2 08 = 2 (/i The one-loop gluon self-energy in the CFL phase has two types of contributions, depending on 

whether the quarks in the loop correspond to singlet or octet excitations, cf. Eq. (23b) of Ref. For the first type, 
both quarks in the loop are octet excitations, and for the second, one is an octet and the other a singlet excitation. 
There is no contribution from singlet-singlet excitations. 

Nonvanishing octet-octet excitations require gluon energies to be larger than 2 08 = 2 0, while octet-singlet excita- 
tions require a larger gluon energy, po ^ 0i + 08 = 3 0. This introduces some additional structure in the imaginary 
parts at po = 30, which can be seen particularly well in Figs.n(d) and (e). 

Some imaginary parts exhibit a peak-like structure at a gluon energy po — E^^ = \/p'^ + (08 + 08)^ = V^O 0, 
followed by a sharp drop for larger energies. The mathematical reason is seen in Eqs. IjAlljl and ljA13|) . where the 
regions po > E^^ and po < E^^ are separated by 8 functions. In the normal phase, the imaginary parts of the gluon 
self-energies actually vanish above po = P- In color-superconducting phases, the imaginary parts do not vanish but 
fall off rapidly. This has already been noted for the 2SC phase 8], and is confirmed here by the results for the CFL 
phase. If there is a peak- like structure at po = E^^, from Eqs. (|A11|) and (jA13|l we expect a similar peak to appear 

at Po — Ep^ = yjp^ 4- (01 -f 08)^ = 5 0. One indeed sees an additional structure at this point, but it is much less 
pronounced since it is located on top of the sharp drop of the first peak. 




In Fig. 121 we show the real parts of the gluon self-energy corresponding to the imaginary parts shown in Fig. ^ 
These are quite similar to the ones in the 2SC phase, cf. Fig. 2 of Ref. |^. An explanation of the various structures 
can be given following arguments similar to those of Refs. |3i In essence, when computing the real part from 
a dispersion integral over the imaginary part, cf. Eq. I)A15|I . a change of gradient in the imaginary part leads to a 
cusp-like structure in the real part. The only function that does not fit this general rule is Re ft'"'. The reason is that 
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this function is computed from the real part of the right-hand side of Eq. Ij4b|l . Note the singularity at po ~ 1.6(f). 
This singularity is caused by a root of the denominator in Eq. Ij4b|l . Pfj^H^" (P) P^, = 0. As shown in Ref. (see 
also Ref. for the 2SC case), this condition determines the dispersion relation of the Goldstone excitations. As one 
expects, for large energies po ^ <j> the real parts of the self-energies approach the corresponding HDL limit. Deviations 
from the HDL limit occur only for gluon energies po ~ 4>. 

We now compute the spectral densities from the real and imaginary parts of the gluon self-energies. When 
Im n'^"(po, p), IniH*(po, p) 7^ 0, the spectral densities are regular and given by 



.00. ^ 1 ImH""(po,p) . 

P (Po,p) = , (5aj 

^ [p2 _Renoo(po,p)]2 + [ImHOO(po,p)]2 



„^ _ 1 ImH*(po,p) 

TT [pI - p2 _ Ren*(po, pW + [IniH*(po, pV^ 



Pbo,p) = Z r^2 ^2 ^.TT..^ , n^^TTf^^ „M2 ■ (5b) 

If ImII°"(po, p) or ImH*(po, p) vanish, the corresponding spectral density has a simple pole given by 



for longitudinal gluons and 



[p2 _ Ren°"(po,p)]po=.oo(p) = (6) 



Jpo=i^*(p) 

for transverse gluons. 



bg-p2-ReH*(po,p)U=^*(p) =0 (7) 




In Fig. 121 we show the spectral densities for longitudinal and transverse gluons in the CFL phase in comparison to the 
HDL limit. Note that there is a delta function-like peak in the transverse spectral density at an e nerg y pq ~ 0.21 rUg. 
This peak corresponds to a collective excitation, the so-called "light plasmon" predicted in Ref. jl2l | (see also [T^). 
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We show the dispersion relation of this collective mode in Fig.0](b). The mass mcoii — 1.35 </> is roughly in agreement 
with the value mcou — 1.362 (f> of Ref. [T^ . As the momentum increases, the energy of the light plasmon excitation 
approaches 2 cj) from below. For momenta larger than ^ 8(j), the location of this excitation branch becomes numerically 
indistinguishable from the continuum in the spectral density above po = 2 0, cf. Fig. 13 Close inspection reveals that 
the dispersion relation of the light plasmon has a minimum at a nonzero value of p ~ 1.33 indicating a van Hove 
singularity. 

In Fig. ^ we also show the dispersion relations for the "regular" longitudinal and transverse excitations, as well as 
for the Nambu-Goldstone excitation defined by the root of Ilfj_,y(P) P" = 8, llj. For our choice of gauge the 
gluon propagator is 4-transverse and this mode does not mix with the longitudinal component of the gauge field Q . 
Therefore, the Nambu-Goldstone mode does not appear as a peak in the longitudinal spectral density, cf. Fig. O We 
finally note that other collective excitations have been investigated in Ref. 




FIG. 4: The dispersion relations for (a) longitudinal and (b) transverse excitations in the CFL phase. The full lines correspond 
to the regular longitudinal and transverse excitations. The dashed lines are for the HDL limit. The dash-dotted line in part 
(a) shows the dispersion relation for the Nambu-Goldstone excitation. The light plasmon dispersion relation is shown by the 
dash-dotted line in part (b). As in Fig. 01 the value of the gap is chosen such that rUg =8 0. 



In conclusion, we have computed the gluon self-energy in the CFL phase as a function of energy and momentum. 
While the imaginary parts of the gluon self-energy could be expressed analytically in terms of elliptic functions (see 
appendix), the real parts had to be computed numerically with the help of dispersion integrals. From the real and 
imaginary parts we constructed the spectral densities. We confirmed the existence of a low-energy collective excitation, 
the so-called "light plasmon" predicted in Ref. [T^. 
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APPENDIX A: CALCULATION OF THE GLUON SELF-ENERGY 

In this appendix, we compute the individual components of the gluon self-energy. According to Eq. (31a) of Ref. 
i. 
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and 
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The projection n°*(P)pi was not explicitly given in Ref. 
steps as for the other components, we find 
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Starting with Eq. (24) of Ref. pf and following similar 
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Here 6 = s-iki — fJ., ki = k and k2 = k — p, where k is the quark three-momentum and p is the gluon three-momentum. 
The octet and singlet gap functions for quasiparticles (e^ = +1) and quasiantiparticles (e^ = —1) are 0i = 0g' and 



0i = 0^', respectively, and the corresponding excitation energies are = v^i+W and = yS,i+4>i- The thermal 

distribution functions are Ni = [exp(ei/T) -I- 1]~^, and Ni = [exp(e,;/T) + 1]^"'^, respectively. In the limit T ^ 0, the 
latter vanish, since e^, > 0. Then the equations simplify to 
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In the following, we evaluate the imaginary part of n"'^ explicitly; the calculation of the other components is similar. 



1. Imaginary parts 

We are interested in the retarded self-energy, so we analytically continue po — *■ Po + i'H in Eqs. (|A2(I . Then, using 
the Dirac identity 
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where V stands for the principal-value prescription, we find 
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Since Eq. ljA4|l is an odd function oipo, Imn(— po,p) = — Inin(po,p)7 the calculation can be restricted to positive 
values of the energy pq > 0. Hence, the first delta-function in each square bracket in Eq. (|A4II can be dropped. 

Moreover, we are interested in gluon energies and momenta po,p fi. Consequently, in order to have a vanishing 
argument of the remaining delta- functions, the quasiparticle energies should not be too large either, e^, fj,. Then, 
only the terms with ei = 62 = -1-1 will contribute; for either ei or 62 = —1, ~ ~ \ki + fJ-\ ^ /i is too far from 
the Fermi surface to make a contribution for po <^ fi. Shifting the integration variable so that ki = k + p/2 and 
k2 = k — p/2, we find using p ^ ^ 
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Choosing p = (0, 0,p) the integration over the polar angle ip becomes trivial. The integral over ^ can be performed 
using the delta-functions. Denoting x — cos 9, the roots of the arguments of the delta- functions in Eq. ljA7|l are (in 
the order of appearance) 
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Then, 
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where y = px/po, $ = (/)/poi = nrin(p/po, yl^^^^VPo ) i ""2 = niin(p/po, v^l — /p^ ) and Wg is the gluon 
mass parameter (squared), to^ = N f g'^ / dn"^ , for A^/ = 3. 



Using the elliptic integrals of first, second, and third kind, 
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and the complete elliptic integrals of the first, K(fc) = F(n/2,k), the second kind, E(fc) = E{TT/2,k) and the third 
kind n(Z,r) — n(7r/2,^,r) in Eq. HA9|I. we obtain the final result 
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where t' = a/ (pp - 9(/)2)/(po ^ 0^) 7 ^ = V'1~^40VPo~' Q^' = arcsin[p/^p^ - 90^ ], a = arcsin[p/(ipo)], ^ = -1 + 

9<jP /p1 and s = 20/po- 

The imaginary parts of the other components of the gluon self-energy can be obtained analogously from Eqs. ljA2b|l . 
"JcJ, and HA2d|l . In addition to Eq. HA5|) we employ 
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l-7e(po-2<^)Je(£;f -po) 



5K(O-(5 + -s2)n0,O 



n=l. 



5i^^(a',i')-(5 + ^s')n(a',;,i') 



In the limit — > 0, we reproduce the HDL hmit 



E(a,i)-iiJl 
Po V 



HmImn°"(P) = Imn!]"(P) , 

Hm Imn*(P) = Imn^(P) , 

linilmn^(F) = Imn^(F) , 

(p—^O 



(A13c) 

(Al4a) 
(A14b) 
(Al4c) 
(A14d) 



2. Real parts 

There are two possibiUties to compute the real parts of the gluon self-energy. Either, one evaluates a principal- value 
integral, cf. Eq. ljA3(l . or one employs the dispersion integral 



1 r 

Ren(po,p) = J 



Imn(a;, p) 



duj '-^^ + C , 

00 -Po 



(A15) 



where C is a (subtraction) constant. If Imn(ti>, p) is an odd function of oj, Imn(— p) = — Imn(aj, p) the dispersion 
integral becomes 



I /I 1 \ 

Re n(po, p) = -V duj ImHoddl^^, p) — + + C 

TT Jq \u;+pq lu-pqJ 

and if it is an even function of Imn(— cj, p) = +Imn(w, p) we have 



I /I 
Ren(po,p) = / dwlmncvcn(w,p) 



1 



(A16) 



(A17) 



where in both cases n(po,p) is assumed to be analytic in the upper complex po plane. 

The values of the constants C^", C*, C^, and are determined by the large-po dependence of the self-energy. 
Thus, it does not matter which color-superconducting phase we consider, and the constants assume the same values 
as for the 2SC phase, = C"* = 0, C* = = m^^ cf. Ref. 0. 
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